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Abstract. Angular momentum plays a crucial role in the formation of 

■ stars and planets. It has long been noticed that parcels of gas in molec- 
• ' ular clouds need to reduce their specific angular momentum by 6 to 7 

, orders of magnitude to participate in the building of a typical star like 

the Sun. Several physical processes on different scales and at different 
^ ' stages of evolution can contribute to this loss of angular momentum. 

. In order to set constraints on these processes and better understand 

C/j ' this transfer of angular momentum, a detailed observational census and 

I characterization of rotation at all stages of evolution and over all scales 

of star forming regions is necessary. This review presents the main 
results obtained in low-mass star forming regions over the past four 
decades in this field of research. It addresses the search and charac- 

■ terization of rotation in molecular clouds, prestellar and protostellar 
^Nj ' cores, circumstellar disks, and jets. Perspectives offered by ALMA are 

, briefly discussed. 

O 

vn 

O ■ 1 Introduction 

cn 

" . [ • Rotation is a ubiquitous phenomenon from the largest to the smallest scales in our 
] Galaxy: the Galaxy itself rotates as a whole, its individual stars spin too - like our 

k>( ■ Sun -, the majority of stars belong to binary systems, and the formation of planets 
; I around stars would not be possible if a certain amount of angular momentum was 
not present at the beginning of the star formation process. The angular momentum 
thus plays a significant role during the process of star and planet formation. The 
typical specific angular momenta measured on different scales and at different 
evolutionary stages, from dense cores in molecular clouds down to the Sun, are 
listed in Table [T] This table shows that a parcel of gas initially located in a dense 
core has to reduce its angular momentum by 6 to 7 orders of magnitude in order 
to participate in the building of a typical star like our Sun. This puzzle has long 
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Table 1. Typical specific angular momenta from dense cores to the Sun 



Object 


J/M 

(c„,2 ,-1) 


References 


Dense cores in molecular clouds 


1021-22 


1 


Protoplanetary disks 


1019-21 


2 


Pre-main-sequence binaries 


1019-20 


3 


Pre-main-sequence stars 


1016-17 


4 


Extrasolar planetary systems (exoplanet(s) 


-t- Star) IQie-is 


5, 6 


Solar system (planets -I- Sun) 




7 


Sun 


1015 


8 



References: 1: Goodman et al. (1993h 2: Williams &: Ciezal (2011'). 3: Chen et al 



(|2007[) . 4: iMathieul (|2004|). 5: lArmstrong et alj (120071) 76: iBerget fc Durrance 
(l2010l) . 7: lAllenI (Il973l) . 8: IPinto et alJ (l201ll) . 



been regarded as the "angular momentum problem" in the field of star formation 
(e.g., ISpitzeij[l97i: iBodenheimeJllQQSl: IMathieul [200I . Obviously, since angular 
momentum is a conserved quantity, this loss of angular momentum has to occur 
by transfer to other particles that will not be incorporated into the star. For 
instance, a fraction of the angular momentum initially present in a dense core can 
be stored by fragmentation into the orbital motion of a binary system. However, 
this fraction is typically on the order of a few percent only (see Table |T]) and this 
process alone cannot solve the angular momentum problem. Magnetic braking 
may also play a significant role in carrying away angular momentum via A lfven 
waves during the early phases of star formation (see, e.g.. iBodcnheimeil ll995l and 
the review of P. Hennebelle in this volume) . 

The specific angular momentum of a dense core has to be reduced by about 
two orders of magnitude during its gravitational collapse until the formation of a 
protoplanetary disk (see Table [T]). During this protostellar phase, the ejection of 
matter in jets and outfiows can partly contribute to this loss of angular momen- 
tum. Later on, the (proto)star-plus-disk system has to reduce its specific angular 
momentum by still a few orders of magnitude to reach the level of our solar sys- 
tem, which appears to have a typical specific angular momentum compared to the 
currently known extrasolar planetary systems (see Tabled]). Finally, the star itself 
has to reduce its angular momentum by one to two orders of magnitude from the 
pre main sequence to the main sequence at the age of the Sun. Star-disk inter- 
actions and magnetized stellar winds are thought to play a significant role in this 
respect (see the review of J. Bouvier in this volume). The solution of the angu- 
lar momentum problem obviously involves different physical processes at various 
stages during the process of star formation. 

In order to better understand when and how these transfers of angular mo- 
mentum occur during the process of star formation and set constraints on the 
physical processes at work, it is essential to collect measurements of rotation on 
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different scales and at different evolutionary stages of the star formation process. 
With this in mind, the aim of this review is to present an overview of the obser- 
vational results that were obtained on rotation in star forming regions over the 
past four decades until February 2013. An earlier review, which also discussed the 
physical proc e sses t hought to be at work on a theoretical basis, was presented by 
Bodenheimer about two decades ago. Since then, many new observational 



results have been collected. These new results, as well as the earlier ones, are re- 
viewed here critically. As a caveat, we stress that this review focuses on low-mass 
star forming regions onlj0. 

The basic tools that can be used to probe rotation in star forming regions are 
presented in Sect. [5] Section [3] addresses the search for rotational signatures on 
large scales in molecular clouds and on smaller scales in dense cores. Sections |4] 
to [6] deal with more evolved stages of star formation: rotation of protostellar 
envelopes is examined in Sect. 21 protoplanetary disks in Sect. [SJ and jets and 
outflows in Sect. |6l The reader will find a short summary of the main results and 
ideas at the end of each section. Overall conclusions are presented in Sect. [7] and 
some perspectives offered by the Atacama Large Millimeter/submillimeter Array 
(ALMA) are briefly described. 

2 How to probe rotation in star forming regions 

This section describes basic tools that can be used to search for signatures of 
rotation in star forming regions. 

2.1 Centroid velocity 

The centroid velocity of a molecular transition observed toward an astronomical 
object is the average velocity of the emitting particles along the line of sight. It 
corresponds to the "systemic" velocity of the object for the line of sight toward 
its center. There are two main traditional ways to derive the centroid velocity. 
For a spectrum with a shape close to a Gaussian, a least-square fitting with the 
3-parameter function 

Tpeake-^'"(2)i^ (2.1) 

yields the centroid velocity Vg, as well as the peak temperature Tpcak and the 
linewidth FWHM. The second method consists in computing the first moment 
of the spectrum: 

N 



i =l 
N 



Vi = '-^ , (2.2) 



^Some aspects o f rotation in high-mass star forming regions are discussed in, e.g., 

iPirogov et all l|2003l l, ICesaroni et al.i (,2007,' ). ,Beltran et al.. t201ll'l . and lLi et aO d2012l'l . 
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Fig. 1. a Projection of the rotation velocity of an object in solid-body rotation along 
the line of sight, for the particular case where the rotation axis is orthogonal to the line 
of sight, viewed from the top. b Configuration when the rotation axis has an inclination 
i to the line of sight, viewed from the side. 



with N the number of channels in the selected velocity range over which the 
emission is detected and Tj the flux density in temperature scale in channel i. If 
the channels are independent, the uncertainty associated to the first moment is 



with (7t the rms noise level of the spectrum. The first method can relatively easily 
handle multiple components along the line of sight but convergence issues are 
encountered when the signal-to-noise ratio becomes too low. The second method 
is mathematically simple and does not evolve any fitting process, but it is more 
difficult to deal with multiple components and the equation may diverge when 

N 

^ Ti is close to zero, for instance in case of a low signal-to-noise ratio. 

i=l 



2.2 How does rotation manifest itself observationally? 

Let's assume an object in solid-body rotation with an angular velocity Cl. Its 



for solid-body rotation, v is proportional to i?, the distance to the rotation axis. 
Let's take two points Pi and P2 along the line of sight, and let's call Pq the point 
closest to the rotation axis along the line of sight and R its distance to the rotation 
axis (see Fig. [T^). The projection onto the line of sight of the velocity difference 




(2.3) 



rotation velocity field is defined by the equation v 



fl X f. fl being uniform 
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between Pi and P2 is 

^d's = I ^ ^ ^ I (t"2 - ri).{v2 ~ Vi) 



r2 - ri 
1 



r2 - ri 



I (^2 - ri).ifl xf2-nxri) 

I r2 ~ ri I V / 
= (2.4) 

The projection onto the Une of sight of the rotation velocity is thus independent 
of the position along the line of sight. This implies that the centroid velocity of 
the spectrum is equal to the rotation velocity of Pq (wccnt — ^R) for a rotation 
axis orthogonal to the line of sight. The proof of Eq. 12.41 is actually valid for 
any inclination i of the rotation axis to the line of sight, and the centroid velocity 
becomes Vcent — ^Rsini for the general case (see Fig.llt). As a result, the velocity 
gradient along the direction perpendicular to the projection onto the plane of the 
sky of the rotation axis is = ^ sin i and is uniform for solid-body rotation. 

If the angular velocity depends on R (but not on the azimuth), the object is 
in differential rotation. For symmetry reasons, the centroid velocity for the line 
of sight with impact parameter R still corresponds to the rotation velocity at Pq 
projected onto the line of sight. Measuring the variations of centroid velocity as a 
function of R thus allows to derive r2(i?) sini. 

A search for rotation signature will start by investigating a map of centroid 
velocity. A solid-body rotation will manifest itself as a uniform gradient of centroid 
velocity, with a magnitude equal to fl sin j. Estimating the mean velocity gradient 
in such a map can thus be done by a simple planar least-squar e fitting with th e 
function t;ccnt = vo+aAa+bA/3, with Aa and A(3 in radian (Goo dman et al.lll993l ). 
The position angle of the direction of the mean velocity gradient is PA — tan~^ (|-) 

and its magnitude is ^'^^+''^ (= risini), with d the distance. The magnitude is 
often given in km 'pc~^ (1 km s~^ pc""'^ = 3.2 10"^'' s~^). 



2.3 Position-velocity diagrams 

Once the direction of the mean velocity gradient is known, it is instructive to an- 
alyze the variations of the centroid velocity in a position- velocity (P-V) diagram 
along this direction. Figure [2^ shows such a diagram for a synthetic protostellar 
envelope in solid-body rotation. The centroid velocity curve is a straight line with 
a slope equal to the velocity gradient. The deviation from a straight line indi- 
cates differential rotation, as illustrated in Fig.[2l3 for an envelope with fl oc r~^-^ . 
The shape of the centroid velocity curve is affected by the finite angular resolu- 
tion of the telescope in the central parts, which prevents drawing any conclusion 
about the velocity field within the inner region of diameter about twice the beam 
width (~ 2HPBW). However, the emission contours in the P-V diagram provide 
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Position (AU) Position (AU) 

Fig. 2. a Position- velocity diagram along the direction of the velocity gradient for a para- 
metric model of a protostellar envelope in solid-body rotation with Q = 4.2 km s^^ pc~^ 
and i = 50° (adapted from iBellochi 120021 ') . The emission (contours) was computed for 
C^^'S 2-1 and a finite angular resolution {HPBW = 3600 AU). The curve marks the 
location of the centroid velocity at each position, b Same as a but for an envelope in 
differential rotation with 57 = 12 (_R/3500 AU)~^'^ km s^^ pc^^. The dashed line shows 
the expected variations of the centroid velocity for an infinite angular resolution. 



additional information about the rotation velocity field even in the unresolved cen- 
tral parts. Fitting a synthetic contour map to the observed one allows to derive 
constraints on this velocity field down to radii smaller than HPBW , provided 
the contribution of other motions (infall, outflow, turbulence) is negligible or well 
known. 



2.4 Angular momentum and rotational energy 

This section lists a number of useful equations for the case of an object in solid- 
body rotation with a power-law density profile p oc r~" and a radius R (see, e.g., 
[Goldsmith &: Arquilla 1985 ). The gradient of centroid velocity is |V?;cGnt| = Hsin? 



(see Sect. 12.21) . For a population of objects with random inclinations, the statistical 
correction to estimate the angular velocities consists in dividing the magnitudes of 
the velocity gradients by (sini) = ^. The moment of inertia is / = |A/i?^ (s^) ' 
the angular momentum J = / 51, and the specific angular momentum j 
jj fi. The rotational energy is i?rot = ^I^"^ and the gravitational energy is: 



j_ 

M 



Esr.v = — ( — ] ■ (2.5) 



R \5-2a 
The ratio of rotational to gravitational energies is 

^ ^rot ^ IR^ fb-2a 
|S,,avl 3 GM 



(2.6) 



Observation of rotation in star forming regions 



7 



For a uniform density (a = 0), we obtain / = |M_R^, j = ^R'^fl, and /3rot = 
i"^ir- ^ centrally-peaked profile with a = 2, / = "^MR^, j = and 
= i . Note that /3i-ot is proportional to — , with pmcan the mean density. 

2.5 Basic procedure to reveal the presence of rotation 

The basic procedure to reveal the presence of rotation requires a tracer being 
optically thin in order to trace all material along each line of sight. Of course, only 
the material at densities (roughly) higher than the critical density of the molecular 
transition is probed. One should in addition keep in mind that the molecule 
abundance may not be uniform due to gas-phase chemical reactions and/or gas 
depletion/release onto/from the dust grain surface. This can affect the derivation 
of the centroid velocity. The advantage of an optically thin tracer is also that its 
lineshape will not be distorted by optical depth effects. 

The first step of the procedure is to search for a one-dimensional, systematic 
pattern in a map of centroid velocity. The second step consists in deriving the 
direction of the mean velocity gradient (via, e.g., planar least-square fitting). The 
final step is to analyze the P-V diagram along the direction of the mean velocity 
gradient. The curve of centroid velocity will indicate if the object is in solid- 
body or differential rotation depending if it is a straight line or not. In addition, 
comparing the emission contours to synthetic P-V diagrams can yield constraints 
on the rotation velocity field at scales even smaller than the beam. 

This procedure is relatively simple but one should keep in mind that, depending 
on the geometry of the source, other types of motion, such as infall, outfiow, or 
shear motions can also produce velocity gradients that could be mistaken as due 
to rotation. 



3 Rotation of molecular clouds and prestellar cores 

This section addresses the search for rotation in the prestellar phase, from the 
scales of molecular clouds down to prestellar cores. 



3.1 Rotation of molecular clouds 

Velocity gradients in tracers such as CO and ^'^CO in emission or H2CO in ab- 
sorption were found in large-scale maps of m olecular clouds i n the 70's. They were 



often interpreted as indications of rotation. iFleck fc ClarkI (|198l[ ) compiled such 
measurements of 13 "rotating" clouds. Their angular velocities $7 normalized to 
the galactic angular velocity J^g ranged from 7 to 300. The early interpretation 
of these high values was that the angular momentum of these clouds was acquired 
from galactic rotation, starting with corotation and then spin-up during cloud con- 
traction. However, this interpretation was questioned by the random orientation 
of rotational axes of early-type stars, binaries, and molecular clouds with respect 
to ^Q. For the cloud sample mentioned above, Vl{R) oc with R the cloud 
radius at which the velocity gradient was measured and p ~ — | (see Fig. 1 of 
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Fleck fc Clarklll98ll) . iFleck fc: ClarkI (|l98ll) interpreted this relation as well as the 



cloud size distribution as resulting from the turbulent properties of the interstellar 
medium, a Kolmogorov cascade with turbulence being sustained on large scales by 
the shearing action of differential galactic rotation. Cloud "rotation" would then 
ori ginate in turbulence vo r ticity . 



Goldsmith fc Arguillal (|l985l) analyzed a sample of 16 clouds without evident 
signs of high-mass star formation. The sample was biased to include only "rotat- 
ing" clouds, i.e. clouds with a clear velocity gradient. The cloud sizes range from 
0.1 to 17 pc with a median of 0.6 pc, and the typical densities are ^ 10'^ cm~^. 
Under the assumption of solid-body rotation, they found specific an gular momenta 



scalin g as R , i.e. fl (x R , consistent with the findings of iFleck fc Clark 



( 198ll ). They interpreted this relation as evidence for loss of angular momentum 



in the process of cloud contraction and fragmentation. They suggested that this 
loss could be due to the redistribution of angular momentum in the orbital motions 
of fragments or to magnetic braking that transfers angular momentum to larger 
scales. Under the assumption of uniform density, they derived /3rot ranging from 
0.04 to 2.5, with a median of 0.25, suggesting that rotation is significant but does 
not dominate the cloud energetics. All the conclusions are based on the assump- 
tion that the velocity gradients do trace rotation. Examining in more details the 
^•^CO 1-0 P-V diagrams of these clouds, the interpretation appears in some cases 
ambiguous: B361 has a P-V diagram suggestiv e of differential rotation alon g the 



axis of its mean velocity gradient (see Fig. 3 of lArguilla fc GoldsmithI 119851 ). but 
the velocity dispersion along the orthogonal axis is as larg e as the velocity varia- 
tions along the axis of the velocity gradient (see Fig. 3b of Arquilla &: GoldsmithI 

Finally, under the assumption that linear velocity gradients observed in their 



sample of 5 giant molecular clouds are due to rotation, llmara fc Blitzj (|201l[ ) find 
that the specific angular momentum of these clouds is smaller than that in the 
surrounding atomic gas out of which these clouds formed. Even more importantly, 
they find that the velocity gradient position angles in the molecular and atomic 
gas are largely divergent, which leads them to suggest that rotation may not be 
the best explanation of the velocity fields observed in giant molecular clouds. 

Overall, all these results suggest that the interpretation of velocity gradients 
on molecular cl oud scales as due to rotation is by no means straightforward and 
robust (but see |Phillipslll999l for an alternative view). 



3.2 Rotation of dense cores 

The presence of rotation on smaller scales correspondi ng to dense cores trace d 



with ammonia (rtHa 10 cm ^) was investigated by iGoodman et al.l (|l993l ). 



Their sample contains 43 cores with radii ranging from 0.1 to 0.5 pc. Significant 
velocity gradients were found in a large fraction of the sample. Some neighboring 
cores have non-parallel velocity gradients (see, e.g., their Fig. 4a). The gradient 
directions do not correlate with the core elongations. The orientations of the 
gradients are random within a complex, for instance for the subsample of 16 cores 
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in Taurus. At least two cores, L1251A and L1251E, have a rotation axis not parallel 
to the magnetic field. Finally, assuming uniform density and not correcting for 
inclination, /3rot is less than 10% for most cores, with a median of 0.03, pointing 
to a minor role of rotation in supporting the dense cores, in agreement with the 
absence of correlation between rotatio n axis and direct i on of elongation. 

The angular momenta derived bv iGoodman et aL (1993) correlate with the 
radii of the NH3 dense cores as R^'^, in a very similar way as for the molecular 
clouds (Sect. 13. 1|) . However, the authors pointed out that, given j oc ilR^, this 
correlation is dominated by R^. The correlation of the angular velocities with 
radius, $7 oc R^'^-^, is actually very weak for their sample. In addition, they did 
not find any correlation between /J^ot and R, which could result from /3rot being 
proportional to ri^/pmoan, ^ not correlating (or only weakly) with i?, and ammonia 
pr obing a small d ynam ic range of densities. 



Caselh et all (j2002h probed rotation of dense cores at even higher densities 



with N2H+ 1-0 (nHa ^ 10^ cm~^). Their sample includes 57 cores, all previously 
mapped in ammonia. On average, the angular velocity derived with N2H+ is 1.6 
times larger than with NH3, but with a large dispersion (1.0). The directions of 
the N2H+ and NH3 gradients are well correlated, and the distributions of P^ot are 
very similar for both tracers ((/3rot) = 0.02 and 0.03, respectively). The planar 
least-square fitting method was also applied on subregions of each core. This 
revealed that many cores have internal variations of magnitude and direction of 
velocity gradient. Simple solid-body rotation thus appears to be rare at the scale 
of dense cores. 



3.3 Probing higher angular resolution 

L183 is a centrally condensed, chemically evolved, high-column density prestellar 
core. Recent interferometric observations in N2H+ 1-0 revealed the presence of a 
velocity gradient on small scales (3300 AU), about 5 times larger than and with 
the same direction as the gradient me asured on larger scales (10300 AU) with 
single-dish telescopes (jKirk et al.ll2009h . The authors interpreted the gradients 
as due to rotation and concluded that the core inner parts span up with rough 
conservation of specific angular momentum. The direction of the velocity gradient 
was claimed to be roughly perpendicular to the projection onto the plane of the sky 
of the magnetic field traced by dust polarization, possibly indicating an important 
role of the magnetic field in the core evolution. However, there are two caveats to 
this result: the velocity gradient is seen only in the northern part of the core, not 
toward the main peak traced with N2H+, and the positions of the continuum and 
N2H+ peaks do not match, suggesting that the N2H+ abundance is not uniform, 
which may affect the measurements of centroid velocities. 

3.4 Complex motions in protoclusters 

The previous sections addressed the search for rotation in relatively isolated dense 
cores. Since many stars form in clusters, it is relevant to investigate the rotational 
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Fig. 3. a Angular velocity as a function of radius for the 5 structures in L1688 with clear 
vel ocity gradient: L1688 as a whole, Oph C-S, E-MM2d, E-MM4, and F-MMl, based 
on iBeUochd l|2002D . b Specific angular momentum as a function of radius for the same 
structures as in a. No correction for inclination was applied and solid-body rotation and 
uniform density were assumed. In each panel, the solid line is the result of a least-square 
linear fitting in logarithmic scales. The slopes are —0.5 and 1.5, respectively. 



properties of prestellar cores located in protoclusters. The protocluster L1688 in 
Ophiuchus con sists of 6 dense clum ps, Oph A to F, within which 60 prestellar cores 
ar e embedded (Motte et al. 19981) . Forty-eight of them were surveyed in N2H+ 1- 
(| Andre et al.ll2007l ). The centroid velocities of the cores trace their motions 



within the protocluster. They reveal a large-scale (^1.1 pc) velocity gradient o f 
magnitude ^1.1 km s~^ pc~^ across the cloud (see Fig. 7 of I Andre et al. 2007 ). 
The position angle (east from north) of the gradient changes from 117° to 186° 
when the cores in Oph B are excluded from the fit. It is thus unclear whether this 
large-scale velocity gradient really traces rotation of the protocluster. The N2H+ 
ce ntroid-velocity ma p of each dense clump reveals its own kinematics (see Fig. 6 
of Andre et al. On the one hand, the velocity structure within Oph A, Bl, 

and B2 is highly complex, with no clear one-dimensional velocity pattern. It is 
thus difficult to conclude anything about rotation in these clumps. On the other 
hand, simple velocity gradients in Oph C and E are consistent with the presence 
of rotation. In Oph F, two unrelated components seem to overlap along the line 
of sight. Each component has a simple velocity gradient, possibly due to rotation. 

The prestellar cores embedded in the clumps with possible signs of rotation 
represent only ~ 25% of the full N2H+ sample of L1688. It is thus difficult to 
draw general conclusions about the rotational properties of prestellar cores in 
protoclusters. However, a closer look at the kinematics in Oph E, wh i ch wa s 
also mapped in ^^^C0+ 1-0, DCO+ 2-1, and DCO+ 3-2 bv I Andre etal] (|2007l ). 



is instructive. The mean velocity gradient over the clump has a magnitude of 
5.2 ± 0.3 km s^^ pc'^ at PA = -108° ± 4°. Fitting the two embedded prestellar 
cores Oph E-MM2d and E-MM4 separately yields velocity gradients of magnitude 
5.2 ± 0.3 km s^^ pc^^ at PA = -126° ± 15° and 7.1 ± 0.4 km s'^ pc'^ at 
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PA = -48° ± 8°, respectively (see Fig. 5.15 of iBellochel |2002| ). The velocity 
gradients around Oph E-MM2d and E-MM4 are significantly not parallel. If they 
are due to rotation, then their misalignment likely results from the turbulent nature 
of rotation in protoclusters. This extends to much smaller s cales the conclu s ion on 
the significant role of turbulence drawn on cloud scales by iFleck fc Clar"3 (|l98ll) 
(see Sect. 13. l| ). This conclusion is further strengthened by the scaling of the angular 
velocity as a function of radius. We find oc i?""-^ and j oc R^'^ for the five 
structures with clear velocity gradient (see Fig. IS]). Even if the statistics is poor, 
this tentative correlation of fl with R in L1688 is close to th e correlation derived on 
larger scales for a sample of clouds bv lFleck fc ClarkI (1981) which was interpreted 
as a signature of interstellar turbulence. 

Assuming solid-body rotation and no correction for inclination, the velocity 
gradients measured for Oph C-S, E-MM2D, E-MM4, and F-MMl imply /3rot rang- 
ing from 2 t o 22% if the d ensity is uniform, or from 0.7 to 7% if p oc (see 
Table 5.7 of iBellochdliooih . Like for the isolated cores mentioned in Sect. 13.21 
rotation is energetically not dominant in the protocluster L1688, down to scales 
of a few 1000 AU. Finally, we note that the rotation period associated with the 
derived angular velocities of the prestellar cores in L1688 is on the order of 10^ yr 
without correction for inclination. This is a f actor 2 to 10 longe r than the life- 
time of the LI 688 prestellar cores estimated bv I Andre et al. (2007). This confirms 
that rotation is dynamically not dominant at the scale of prestellar cores in this 
protocluster. 



3.5 Peculiar motions in starless cores 

This section addresses two cases further illustrating the difficulties in probing ro- 
tation in starless cores. L1506C is a large, low-density, starless core embedded in 
a filament in Taurus. Its inner parts (r < 0.15 pc) are characterized by a very low 
level of turbulence (crturb < 47 m s~^) and a high level of C^^O depletion, unex- 
pected given the low density of the core (IPagani et al.ll2010l ). The authors proposed 
that this high level of depletion is related to the low level of turbulence promot- 
ing dust coagulation, which in turn may decrease the desorption efficie ncy. An 
even m ore puzzling property of this core is its kinematical structure: .Pagani et al.l 
( 2010l ) report a velocity gradient in ^^CO with a direction opposite to the velocity 



gradients traced by C^^O and N2H"'". Their detailed radiative transfer analysis 
shows that the inner parts are contracting (v-mf = 0.11 km s~^) and rotating with 
fl oc r^^'^, and the o uter parts are expan ding {vcxp — 0.09 km s^^) and rotating 
in opposite direction. iPagani et al.l ([20i3) interprete these properties as evidence 



for decoupling of the inner core from the external parts, and suggest that L1506C 
is a prestellar core in the making. The peculiar velocity fields (infall/expansion, 
counter-rotation) remain however to be explained physically: the authors suggest 
that they may result from oscillations produced by magnetic torques. Alterna- 
tively, the velocity gradients may not trace rotation. 

The prototypical, isolated, star le ss core B68 is sup ported by thermal pressure 
( Alves et al...2001:.Lada et al.ll2003[ ). Lada et al. (2003) reported an approximately 
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east- west velocity gradient, which they interpreted as rotation with ffmt ^ 4% 
Howev er, the presence of a "bullet" to the south-east of B68 led lBurkert fc Alves 
(|2009l) •;o propose that B68 is undergoing a collision with this small core, which 
they qualitatively repr oduced with hydrodyn amic simulations. This scenario was 
recently followed- up bv lNielbock et al" ( 2012 ) who analyze larger maps of centroid 
velocity. They establish a connection between the systemic velocity of the putative 
colliding small core and the velocity gradient seen across B68 (see their Fig. 17). 
On even larger scales, their ^^CO 2-1 centroid- velocity map shows that the velocity 
structure of B68 is related to a larger-scale gradient along a large-scale underlying 
filamentary structure (see their Fig. 18). This velocity gradient could be due to 
streaming motions of cores along the filamentary structure rather than rotation. 
This would support the collision interpretation for B68. 



3.6 Reliability of angular momenta derived observationally 

The rotation measurements reported in the previous sections are based on two- 
dimensional centroid-velocity maps. The loss of the third spatial dimension inher- 
ent in astronomical observations may bias the interpreta tion of vel o city g radients. 



ent m astronomical observations may bias the interpreta tion ot vel o city g radients 
and in particular the calculation of angular momenta. iDib et all tOldi investi 



gate this issue by means of three-dimensional magnetohydrodynamic simulations 
of isothermal, self-gravitating clouds with decaying turbulence. The rms Mach 
number of turbulence in their simulations is similar to those of the Ophiuchus 
(L1688) and Perseus molecular clouds. The dense cores generated in these sim- 
ulations show a variety of morphologies, from roundish to filamentary (see their 
Fig. 3). The authors compute projected maps of centroid velocity to compare 
the outcome of the simulations to observations. As their Fig. 9 shows, turbulence 
generates ve locity gradients qualitatively similar to those measured in molecular 
cloud cores. Dib et ahl ( 2010( ) compute the distribution of 3D specific angular mo- 



menta (j3d) of their synthetic magnetized cores. They find a median about 5 to 10 
times lower than the median of angular momenta (j2D,obs) derived observationally 
from two-dimensional centroid-velocity maps of NH3 and N2H+ dense cores (see 

To understand this discrepancy, they also compute the specific angular mo- 
menta of the synthetic cores by following the observational procedure (j2d), i-e. 
measuring the mean velocity gradient in projected maps of centroid velocity and 
assuming solid-body rotation and uniform density to compute the angular mo- 
mentum. They find ^20 one order of magnitude larger than jan for their synthetic 
cores, i.e. ^20 values consistent with the observed ones (j2D,obs)- The "rotational" 
properties of the synthetic cores are thus very similar to the observed ones, but 
this analysis suggests that the assumption of solid-body rotation and uniform 
density to compute the specific angular momentum of a prestellar core based on a 
two-dimensional map overestimates the true value of the angular momentum. As- 
suming centr ally peaked densit y profiles is not sufficient to remove the discrepancy. 
Interestingly, Dib et ahl ( 2010[ ) find a correlation between jsd and i?3D somewhat 



steeper than between J2d and i?2D (exponent 1.8-2 vs. 1.1-1.3). 



Observation of rotation in star forming regions 



13 



A puzzling result of this work, however, is that the synthetic and observed 

cores have similar distributions of /Spot (/^rot.sD for the former and /3rot,2D,obs for 

the latter) while Prot,3B < /5rot,2D,obs could have been expected a prio ri given that 

. irnt.m < Jrot,2D.obs- A key result to solve this puzzle is provided bv lOffner et al.l 

( 2008f) who did a similar analysis for a population of non- magnetized synthetic 

cores. They obtain 4^ ~ 0.01-0.1 for the synthetic cores (see their Fig. 5), in 

rough agreement with the findings of lDib et al. I (l2010l) . In addition, they compute 

for the synthetic cores and obtain a median value close to 1, which explains 
''¥^ 1 I 1 

whv lDib et al.l (|201C1I) find a good agreement between /3rot,3D and /3rot,2D,obs- The 
reason why o"*'^" ~ 1 while 4^ ~ 0.01-0.1 is unclear. We speculate that it is due 
to the scalar nature of /3rot,3D as opposed to the vectorial nature of the angular 
momentum. Two parcels of gas with opposite J will not contribute to the total 
angular momentum while they will both contribute to the "rotational" energy. 
This illustrates that /3rot certainly includes more than pure rotational energy. 



3. 7 Summary 



Given their scaling properties, the velocity gradients measured in molecular clouds 
and interpreted as rotation likely originate in turbulence vorticity. At the scale 
of dense cores, the rotational energy is only a few percent of the gravitational 
energy. Rotation is thus energetically not dominant. The direction of rotation 
axis is not correlated with the core elongation or the magnetic-field direction. The 
correlation between the specific angular momentum and the radius {j oc i?^'^) 
is similar to the one found on larger scales for molecular clouds. Interpreted in 
terms of rotation, it suggests a loss of angular momentum during the contraction 
process of dense cores. This loss could be due to magnetic braking, gravitational 
torques, or result from the transfer of angular momentum into the orbital motion 
of fragments. However, as for clouds, this behaviour may simply result from the 
properties of interstellar turbulence. The interpretation of velocity gradients in 
clouds and dense cores as rotation is therefore not unique. In any case, although 
often used as a zeroth-order approximation, solid-body rotation appears to be rare, 
which, in combination with the two-dimensional nature of astronomical data, likely 
results in overestimating the specific angular momenta and the ratios of rotational 
to gravitational energies. Finally, the velocity structure in protoclusters or "non- 
isolated" cores is often complex, resulting from, e.g., confusion along the line 
of sight, shear motions, or streaming motions along filamentary structures. In 
protoclusters, turbulence seems again to play a major role in the production of 
velocity gradients. Probing rotation at the prestellar stage is thus often more 
challenging than one could naively expect. 
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4 Rotation of protostellar envelopes 

4.1 General properties of protostars 

Protostars are systems in the main accretion phase, with a stellar embryo - i.e. 
an hydrostatic object - at their center, accumulating mass from a collapsing enve- 
lope and/or accreting circumstellar disk. This accretion process is systematically 
accompanied by ejection of matter in the form of bipolar jets and outflows. Two 
types of protostars have been defined observationally: Class protostars have 
most of their mass still in the envelope ( Mp^y > M^,) w hile Class I protostars have 
accreted most of their mass (A/o„v < M^) (|Andre et al.|[r993.. ,2000 '). In the Class 
phase, the envelope is still prominent and should retain memory of its initial con- 
ditions, and in particular have similar properties as in the prestellar phase. In 
the Class I phase, the system is dominated by a star-disk system, with a residual 
envelope. 

Studying protostars is highly relevant for understanding star formation, and 
in particular the evolution of angular momentum, because a protostar will form 
a star. This is not necessarily the case fo r starless cores, som e of them being 



maybe only transient structures (see, e.g., iBelloche et al.ll201ll for a discussion 
of the fate of starless cores). The advantage of protostars over starless cores for 
searching for rotation is that the position of the system "center" is known. In 
addition, jets and outflows, wh ich arc thou ght to be driven by magnetocentrifugal 



acceleration (e.g.. lK6nigl fc Pudritz 2000: Shu et al.ll2000 ). define a "natural" axis 



for rota tion a priori. In addition , they can give clues about the inclination of the 
system (jCabrit fc BertoutI [l990[) . which is needed to derive the angular velocity 



and angular momentum. 



4.2 Rotation and outflow axis 



As stated in Sect. 14. ll we naively expect to find velocity gradients tracing rotation 
in protostellar envelopes in the dir ection perpendicular to the outfiow axis. In 
practice, this is not often the case. ICurtis fc Richeii (|201l[ ) surveyed starless and 
protostellar cores in Perseus in C^^O 3-2. Seven Class cores have a simple bipolar 
CO outflow and a significant C^^O velocity gradient. For two of them, the velocity 
gradient is roughly orthogonal to the outflow axis, as we expect for rotation, but 
for 5 of them, it is roughly parallel. This may indicate that the rotation axis 
changes direction from large scales (nHa ^ 10^-10^ cm~^) to sma ll scales where 



the je t /outflow is launched (at most a few AU, see Sect. |6] and, e.g.. lFerreira et al 
20061 ). Alternatively, C-'^^O 3-2 observed with single-dish telescopes may be a 
poor probe of envelope rotation. C^^O is indeed known to suffer from depletion at 
high density. Entrainment by the outfiow may in addition contaminate the P-V 
diagrams. In high-density protoclusters, C^^O 3-2 may also be too sensitive to 
the ambient cloud. 



Tobin et al.l ()2011l) surveyed 17 isolated, mostly Class but also a few Class I, 



protostellar envelopes in N2II+ 1-0 and NH3 (1,1) with single-dish telescopes and 
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interferometers. The single-dish maps reveal 11 sources (out of 16) with a velocity- 
gradient direction lying within 45° of the normal to the outflow, and 12 sources out 
of 14 observed interferometrically have the same property (see their Fig. 27). In 
addition, the directions of the single-dish and interferometric gradients are found 
to be generally consistent (see their Fig. 28). The average magnitude of the single- 
dish velocity gradients is 2.2 km s~^ pc^^, and 8.6 for the interferometric data. 
From these results, it is tempting to interprete the velocity gradients as due to 
rotation, with spinning-up from large (~ 10000 AU) to smaller 1000 AU) scales. 
However, there is a significant spread of position angle relative to the normal to 
the outflow, which may indicate that the velocity gradients do not trace pure 
rotation. But if they do, then, again, these offsets imply that the orientation of the 
rotation axis changes from scales of ^ 10000 AU to scales of < a few AU where the 
jets/outflows are launched. The authors suggest also that the rotation signature 
may be affected by the a symmetry of the en velopes, most of them being indeed 



fo und to be asymmet ric (jTobin et al.ll2010al ). Note also that both tracers used 



bv lTobin et al.l ( 201 if ) do not peak on the central protostar in general, suggesting 
th at they do not trac e well the innermost parts (r < 1000 AU). 

Chen et al.l pOOTI ) carried out a similar interferometric survey of 8 Class 
pr otostars in NoH"*" 1-0. Only two of these sources, which are also in the sample 
of iTobin et al.l (j201l[ ). have a velocity gradient close to the normal to the outflow. 
Combining both samples, only 12 out of 19 Class envelopes have this property 
on scales ~ 2000-8000 AU. We conclude from this that even N2H"'" 1-0 is not 
an ideal probe of rotation in protostellar envelopes, or, again, that the angular 
momentum on scales at which the jet/outflow is launched (< a few AU) has a 
significantly different direction than on scales of a few thousand AU. In the latter 
case, it would be tempting to conclude that turbulence vorticity, rather than well- 
ordered rotation, still dominates the velocity field on scales of a few thousand AU. 
It would then be interesting to check whether the protostellar envelopes with a 
velocity gradient not perpendicular to the outflow axis are also those with a higher 
level of turbulence. 

Finally, in this context, it is worth mentioning that the TADPOL polarization 
survey performed toward a sample of 27 protostellar sources with CARMA found 
a random orientation of the magnetic fields with respect to the outflow axes on 



scales of ~ 1000 AU ( Hull et al.ll2012f ). This and the results mentioned in the 



previous paragraphs suggest that the naive picture of the rotation, magnetic field, 
and outflow axes being aligned in a protostellar envelope is rare at this scale. 
For completeness, it would be interesting to investigate if there is a degree of 
correlation between the magnetic field direction and the direction of the velocity 
gradients thought to trace rotation in these protostellar envelopes. 



4.3 The test case I RAM 04191 



IRAM 04191-1-1522 - hereafter IRAM 04191 - is a very young Class protostar 
(1-3 X lO'' yr) located in Ta urus. It has a pro minent envelope 1.5 Mq) and 



drives a collimated outflow ( Andre et al.lll999j ). Its internal luminosity is very 
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Fig. 4. Position-velocity diagram of C'^'^S 2-1 along the direction perpendicular to the 
outflow axis in TRAM 04191. a Synthetic contour map of the "best-fit" collapse model 
with differential rotation, b Observed contour map. The black dots with error bars are 
the measured centroid velocities. The red curve in each panel is the synthetic centroid- 
velocity curve. The angular resolution is 25.5" (HPBW), i.e. 3600 AU (figure taken 
from iBelloche et aLll2002h . 



low (0.08 L0. iDunham et al-lliooi ). IRAM 04191 belongs to the "well-behaved" 
protostellar envelopes in terms of velocity gradient: it harbors a regular, large- 
scale velocity gradient the direction of which is paral lel to the envelope m ajor 



axis and orthogonal to the outflow axis (see Fig. 2 of IBelloche et al.ll2002l ). Its 



centroid-velocity curve along the direction of the mean velocity gradient is not a 
straight line, but it is centrosymmetric, which is a strong evidence for rotation 
(see Fig. \^). The magnitude of the velocity gradient is 7-10 km s~^ pc~^ at 
r - 2800 AU 1.3 km s~^ pc~^ a t r ~ 7000 AU, and < 0.8 km s~i pc^^ at 



11000 AU (jBelloche et al.ll2002l) . The envelope is thus in differential rotation, 



the inner parts rotating faster than the outer parts. The observed P-V diagram 
is well reproduced by a parametric model of a differentially rotating, infalling, 
spherical envelope, with CS (and C'^^S) depleted toward the center (see Fig. H^). 
The model fits well several transitions of CS and C'^^S and strong constraints 
o n the infall and rota tion velocity fields could be derived (see Figs. 14 and 12 
of IBelloche et ah 2002). These kinematica l features are in qualitative agreement 
with a model of Basu fc MouschoviasI ( 1995 ) of a supercritical magnetized envelope 



collapsing and detaching from its subcritical environment (see IBelloche et al.ll2002 
and Fig. 0. 

The differential pattern seen in the P-V diagram of IRAM 04191 is not unique: 
the P-V diagram of CB230 shows a similar shape for instance (see Fig. 24 of 
Tobin et al.ll201ll ). B ut only 2-3 protoste llar envelopes have a similar P-V diagram 



out of 17 observed bv lTobin et al. (2011). Therefore the conclusions obtained for 
IRAM 04191 cannot be generalized to all Class protostars. 
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Fig. 5. Comparison of the profiles of angular velocity (a) and radial velocity (b) of model 
8 of lBasu fc Mouschovi a's (1995) with the profiles (c and d) derived for IRAM 04191 (see 
caption of Fig. 4-9 of lBellocha.2002i for more details). 



On smaller scales, the angular velocity of the envelope of IRAM 04191 keeps 
increasing. BcUochc & Andre (2004a) derive a velocity gradient of 26 km s~^ pc~^ 
at r 2000 AU from N2H+ 1-0 interferometric observations. However, the N2H+ 
emission map presents a hole toward the center and the column density profile 
is consistent with N2H"'" being absent from the gas phase for r < 1600 AU (see 
Figs. 1 and 4b of iBelloche fc Andrei l2004bl) . Because there is no evidence for 
strong C^^O desorption. iBelloche fc Andrl (|2004b[ ) conclude that N2H+ suffers 
from depletion at densities higher than 5 x 10^ cm~'^. Based on in dependent 



N2H+ 1-0 interferometric measurements with similar angular resolution, iLee et al 
( 2005( 1 report a linear velocity gradient for r < 1400 AU. They conclude that the 
innermost parts of the envelope are in solid-body rotation, indicating that the 
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Fig. 6. Position velocity diagrams of a spherical or axisymmetric envelope (top row) 
and a filamentary envelope (bottom row), with rotation only (right column), infall only 
(middle column), and rotation and infall (left column) (from Flbbin et al ] |2012al ). 



infalling material loses angular momentum, possibly due to magnetic braking like 
in the collapse simulations of magnetized singular isothermal toroids bv lAUen et al 



However, an alternative (more likely?) explanation is that the depletion 
of N2H"'" from the gas phase significantly affects the P-V diagram. If there is no 
molecule emitting below a radius tq, then the P-V diagram for R < vq will trace 
the velocity field of the shell of radius tq only, which, in projection, will appear 
as a linear velocity gradient. A tracer less sensitive to depletion is thus necessary 
to unambiguously probe the rotational properties of th e innermost parts of t he 
IRAM 04191 envelope. H2D+ may be a good candidate (jWalmslev et al.ll2004l ). 



4.4 Do velocity gradients really trace rotation of protostellar envelopes? 

The protostellar envelopes studied by iTobin et al.l (|2010a[ ) are often flattened . 
They even find indications that some of them are filamentary. iTobin et al.l (|2012al ) 
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Fig. 7. Specific angular momentum as a function of radius for a sample of sources in 
Taurus. The squares are NH3 dense cores ([Goodman et al.lll993h . The stars are mea- 
surements in TRAM 04191 at four different scales (|Belloche et al ] |2002l : lBelloche fc Andrei 
l2004al ). The filled circles are rotating, infalling envelopes and the triangles are "rotation- 
ally supported disks" , all from lOhashi et al] ([l993) (adapted from lOhashi et aTlllQQTl and 
iBelloche et al] I2OO2I ') . 



argue that infall in such prolate structures at intermediate inclination would pro- 
duce a gradient in centroid-velocity maps that can mimic (differential) rotation 
(see Fig. [p. As a proo f of c oncept, they use an analytic model of rotating coUapse 
based on lUlrichI ( 197d ) and Cassen fc MoosmanI (1981), where the self-gravity of 
infalling gas is neglected. The particles follow ballistic trajectories around a cen- 
tral gravitating mass. The model is truncated t o mim ic a filamentary structure. 
This model is applied to 5 sources of Tobin et al.l ( 2011 ) with a well-ordered veloc- 
ity field. Most P-V diagrams can be reasonably well fitted with such a model of 
collapsing filament. The agreement between model and observations is not a proof 
that the kinematics is effectively dominated by infall in a filamentary structure, 
but it shows that this idea is a viable concept. If envelopes are really filamentary, 
then the velocity gradients could probe infall rather than rotation. In this case, 
it would become necessary to go to smaller scales to trace rotation and derive 
meaningful constraints on the distribution of angular momentum. 



4.5 Specific angular momentum in protostellar envelopes 

In this Section , we assume t hat velocity gradients do trace rotation in protostel- 
lar envelopes, lohashi et al.l (|l997i ) compiled measurements of velocity gradients 
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in a few protostellar envelopes and circumstellar disks of the Taurus molecular 
cloud made in the 90's. The specific ang ular momenta of these systems are com- 
pared to those of the NH3 dense cores of iGoodman et alj (|1993l ) in Fig. [71 which 
can be interpreted as an evolutionary diagram. The measurements obtained for 
IRAM 04191 (see Sect. 14.3) ) are also displayed. The local specific angular momen- 
tum is approximately constant for infalling envelopes and circumstellar disks in 
Taurus (about 10"'^ km pc or 3 x 10^° cm^ s~^, see dashed line in Fig. [7]), 
while the angular velocity is approximately constant for the dense core regime (see 
solid line in Fig. [7]). The transition between the dense core j-R correlation to a 
regime of constant specific ang ular momentum occu rs at a radius of ^ 5000 AU 
or ~ 0.03 pc. As pointed out bv lOhashi et al. (1997), this appears to be the scale 
below which the collapse is dynamical with conservation of angular momentum. 
The interpretation of the inner envelope of IRAM 04191 being a magnetically su- 
percritical core detaching from its subcritical environment fSect. l4?3l) fits relatively 
well into this patter A more complete version of this diagram will be presented 
and discussed in Sect. 17.21 



4.6 Summary 

The velocity structure of protostellar envelopes is often complex. In about 50% 
of the cases, the velocity gradient in the envelope is not orthogonal to the outflow 
axis. This suggests that either velocity gradients do not reliably trace rotation (on 
scales 1000-10000 AU), or the direction of the rotation axis changes from large 
(~ 10000 AU) to small (< a few AU) scales where the jets/outflows are launched. 
Given that protostellar envelopes are often asymmetric, the velocity gradients 
may also trace infall if the envelopes have a filamentary structure. Therefore, 
velocity gradients do not unambiguously trace rotation in protostellar envelopes. 
However, if the rotation signature in P-V diagrams can be trusted, there are good 
indications for differential rotation, with spinning-up of the inner parts, to occur 
below R ^ 0.05 pc (e.g., IRAM 04191). In addition, there seems to be in Taurus 
a characteristic scale of i? ~ 0.03 pc below which the specific angular momentum 
is conserved, probably as a result of dynamical collapse. 



5 Rotation of protoplanetary disks 

5.1 General properties of circumstellar disks 

The presence of disks aroung young pre-main-sequence (PMS) stars was initially 
inferred from the infrared excess emission in the spectral energy distribution of 
these stars. The flattened morphology of these putative disks was later conflrmed 



■^Beware, however, that the datapoints for IRAM 04191 in Fig.[7|represent the profile of specific 
angular momentum for an envelope at a given time, while the other points rather correspond 
to the evolution in time of the specific angular momentum of protostellar systems, which is 
conceptually different. A direct comparison is therefore not straightforward. 
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by imaging, in part icular as shadows in front of a bright background in the optical 
(see, e.g.. Figs. 1 of lMcCaughrean et al.lll998l and lSmith et al.ll2005l) . The lifetime 
of disks a round solar-type PMS s tars is on the order of 1-10 Myr with a median 
of 3 Myr ( Williams fc Cieza 20 111). Their mass ranges from 0.5 to 50 Mj^p with a 
median of 5 Mjup = 0.005 Ma, (I Williams fc C iezall201l[). Given the high detection 

lUfc , it is 1 



rate of exoplanets around stars (e.g.. lMavor et al. 20111). it is reasonable to consider 



the circumstellar disks around PMS stars as protoplanetary. 



In the Class II phase, the median ratio of disk to star masses is 0.9% ([Williams fc Cieza 



20111 ). Therefore the disk self-gravity is negligible and we expect their kinematics 
to be dominated by Keplerian motion, i.e. Vrot = \J ■ Centroid-velocity maps 
of such disks are expected to have a very characteristic bipolar morphology (see. 



e.g.. Fig. 2 of lGuilloteau et al.ll2006[) . Since Keplerian disks exist in the Class II 



phase (see Sect. 15. 5p . they must have formed during the Class O/I protostellar 
phase. The next three sections will try to address the questions of when a disk 
appears in the process of star formation and when Keplerian rotation sets in. 



5.2 Disk formation 



The centrifugal radius is defined as the radius where gravity is balanced by the 

2 

centrifugal force For a collapsing, non- magnetized, singular isothermal sphere 
initially in solid-body rotation, the centrifugal radius is found to grow with time 
as: 

r-i \2/ \-8/ \3 

(5.1) 



i?c(AU) - 39 



10-14 rad s" 



0.2 km s" 



1 Mq 



with m*+d — 0-975^t and a the sound speed (jTerebev et al.lll984l ). Under these 
restrictive and idealized conditions, rotationally supported disks are expected to 
form from the beginning of the Class phase on with a radius growing as t^. 
As an illustrative application, let's compute the expected centrifugal radius of 
IRAM 04191 in the framework of this model. We use the velocity gradient mea- 
sured at 2000 AU to estimate ft (see Sect. 14. 3| ). For a temperature of 10 K and an 
age of 1-3 X lO'' yr, we find Rc = 3-80 AU. The constraints derived on the kinemat- 
ics of the envelope imply an up per limit to the centrifugal radius of 400 AU (see 
Fig. 15 of iBelloche et al.l 120021 ). In addition, the weak interferometric detection 
in continuum emission at 1.3 mm implies either a disk mass lower than 0.001 Mq 
or a disk r adius smaller than 1 AU if the emission is optically thin or thick, re- 
spectively (jBelloche et al.ll2002l ). Depending on the exact age of IR AM 04191, the 
latter upper limit is marginally consistent with the prediction of iTerebev et al 



The non-magnetized picture of disk growth described in the previous paragraph 
has been challenged in the past decade by the advent of 3D MHD simulations. In 
ideal MHD, for magnetization levels comparable to those observed in dense cores, 
these simulations predict a "magnetic braking catastrophe" which prevents the 
formation of a disk when the magnetic field is parallel to the rotation axis (e.g.. 
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Mellon fc Lil 2008 : Hennebelle fc Fromane 2008t see also the axisymmetric simu- 
lations of lAUen et al.l 120031 ) . When the initial magnetic and rotation axes are not 
aligne d, disk formation can occur depending on the initial mass-to-magnetic-fiux 
ratio (jHennebelle fc Ciardill2009l: I Joos et al.ll2012l) . Non-ideal MHD effects such as 
ambipolar diffusion, Ohmic dis sipation, an d the Hall effect, have also been investi- 
gated recently (e.g IShu et al. 2006: Mell on fc Lill2009l: iKrasnopolskv et al IbOlOt 
iDapp fc Basull2O10HMachida et al.ll201ll:lKrasnopolskv et al.lbOllHLi et al.ll201l[) . 



Ambipolar diffusion does not change the disk outcome, but Ohmic dissipation and 
the Hall effect may (or not) enable the formation of rotationally supported disks 
in the Class p hase, albeit maybe small ones only (< a few AU, see discussion 



in iLi et al.ll201ir) . This is currently an active field of research and the conclusions 
are not well settled. Observational constraints on the formation of disks in the 
Class phase are highly needed. 



5.3 Class protostellar disks (?) 

Several interferometric surveys were performed in dust continuum emission to 
search for disks around Class and I protostars. The PROSAC survey was per- 
formed at 1.1 mm with the SMA with an angular resolution of 1-3" (i.e. not 
sufficient to resolve most disks), with complementary single-dish maps obtained at 
870 /im with SCUBA on JCMT. The sample i ncludes 10 Class and 9 Class I pro- 
tostars at distances ranging from 125 to 325 pc (jj0rgensen et al Ilooi). The respec- 
tive contributions of the envelope and the (unresolved) "disk" to the continuum 
emission was estimated using a parametric model of spherically-symmetric proto- 
stellar envelopes with a power-law density structure. Assuming 3 K for the disk 
emission and the same dust opacity coefficient as for the envelope. I Jorgensen et al ' 



( 2009f) derived a median disk mass of 0.09 Mq for Class and 0.011 Mq for Class I 



objects. After corrections for non-uniform temperature, they obtain the same me- 
dian disk mass for Class and I protostars (~ 0.04 Mq). Taking grain growth into 
account would lower the Class I median mass with respect to the Class one. It 
appears from this study that (massive) disks are formed early in the Class phase 
and that the infalling material from the envelope must be rapidly transported 
thr ough the disk onto the star. 



Enoch et al. ( 201l[ ) surveyed 9 Class and 3 Class I protostars in Serpens 



with CARMA at 1.3 mm with an angular resolution of 1". They estimated 
the "disk" masses from the long-baseline (50 kA) flux, arguing that the envelope 
emission contributes to at most 30% and can be neglected. They detected 6 Class 
sources and obtained a median disk mass of 0.15 Mq. Based on this small sample, 
they did not see any obvious systematic variation of the d i sk ma ss from the Class 
to the Class I phase, i n agreement with Jorgensen et al. ( 20091) . 



Maurv et al. I (l2010l) compare the results of their sub-arcsecond resolution {HPBW = 0.3-0. 



dust continuum survey of 5 Class protostars in Taurus and Perseus performed 
with the IRAM Plateau de Bure interferometer (PdBI) to the outcome of hydro- 
dynamic simulations of a massive (0.7 M©) disk around a 0.7 M© sta r, which frag- 

ments rapidly and may thus exist only during the Class phase (|Stamatellos fc WhitworthI 
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20091) . They conclude that the data are not consistent with the model and that 



none of their sources can harbor such a massive disk. However, they do not exclude 

the presence of 0.1 Mq disks^ 

Although the results of |j0rgensen et al. (|2009t) and lEnoch et all (|201ll ) tend 
to suggest that ~ 0.1 Mq "disks" are common in the Class phase, a significant 
number of Class protostars do not show any e vidence for the pre sence of such 
disks: Md isk < 0.02 Mq for 3 so urces in Serpens dSnoch et al.ll201ll). < 0.0 1 Mq 



for L483 (jjorgensen et all 1200 9'). < 0. 006 Mq for L 72 3 (jCirart et all bool) . and 



< 0.001 Mq for IRAM 04191 (BcUoc he et al.1 120021) . L1521F is even weaker than 



fM 



IRAM 04191 in the survey of iMaurv et al.l (|2010[ ). suggesting that its disk, if any, 

is even less massive. 

Although |j0rgensen et al" ( 2009[ ) and Enoch et al. ( 2011 ) claim detecting ~ 
0.1 Mq "disks" in their samples of Class sources, it is important to keep in 
mind that the emission was not resolved by these observations and that strong 
assumptions about the structure of the envelope - spherically-symmetric envelope 
with power-law structure for the former study, no envelope contribution at 50 kA 
for the latter - had to be made in order to derive the excess emission attributed 
to a "disk" . IChiang et al . (2008) compared the emission of 4 Class sources 
observed by iLoonev et al. [ (l2000l) at sub-arcsecond resolution {HPBW ^ 0.6") 
with BIMA - 3 being also in the PROSAC sample - with a model of a collapsing, 
magnetized envelope ( Tassis k. Mouschoviad 20051 ) . They obtain good fits for all 
their sources without the need for an additional circumstellar- disk component 
at the 90% confidence level, which contradicts the conclusions of J0rgensen et al.l 
(200i) for the sources that both studies have in common. This shows that knowing 
the inner structure of the env elope, which departs from a power-law in the model of 
Tassis fc Mouschoviaj (2005), is crucia l to derive reliable constraints on the mass 



of a putative disk. Chiang et al. ( 20081 ) mention that models with a disk mass up 
to ^ 0.1 A/0 in addition to the magnetized envelope give acceptable fits a ls o. 

The unresolv e d com ponents interpreted as disks by IJorgensen et al. I llooi) 
and Enoch et al.l (2011) may actually be infalling, magnetized pseudo-disks (e.g.. 



Hennebelle fc Fromanal2008[ ) , i.e. structures that become flattened because of the 
magnetic field rather than rotation. The observed unresolved components coul d 
also be related to outflow cavities rather than circumstellar disks ( Joos et al]|2012 ). 
Higher-angular resolution is necessary to better resolve the continuum emission 
and allow for a detailed comparison to predictions of collapse models. In addition, 
probing the velocity field of these compact structures will be crucial to distinguish 
between the alternatives mentioned above: for instance, a magnetic pseudo-disk 
should be infalling while the velocity field of a rotationally supported disk should 
be dominated by rotation. 

One of the highest-resolution study of the kinema t ics of the inner parts of a 
Class protostellar envelope was done by [Ch oi et al.l (|2010t ). They mapped the 
Class binary system NGC 1333 IRAS 4A in NH3 (2,2) and (3,3) with VLA with 
an angular resolution of 0.3", i.e. 70 AU at the distance of 235 pc which they 
assume for Perseus. They detect a compact emission of deconvolved size 130 x 70 
AU around the component 4A2 and find a clear velocity gradient perpendicular 
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to the outflow axis (see their Fig. 2). They use a parametric model of a thin disk 
assuming a power-law distribution for the NH3 line brightness to fit the velocity 
field of this compact structure. They reproduce well the position-velocity diagram 
taken along the direction perpendicular to the outflow assuming a rotation velocity 
Vrot oc r~°'^. Assuming the structure is in Keplerian rotation, they derive the mass 
of the central star (0.08±0.02 Mq) and then its age (5 x 10"* yr) from an estimate of 
the mass accretion rate (2 x 10~^ Mq yr~^) based on the measured accretion lumi- 
nosity and assuming R^, = 2 Rq . This age is somewhat longer than the dynamical 
time associated with the outflow (1.5 x 10** yr, Knee & Sandcll 2000), but is not 
inconsistent given the large uncertainties on the luminosity and the stellar radius. 
This could be the first reliable observational evidence for a rotationally supported 
disk in the Class phase. We note however that the mass infall rate measured in 
the collapsing envelope is about 50 times higher than the ma ss accretion rate men- 
tioned above ( Di Francesco et al. 2001; Belloche et al. The mass accretion 



rate may thus have been underestimated and radial motions in the "disk" may 
still be significant, implying in turn that the measured rotation velocities may be 
sub-Keplerian and not Keplerian. In any case, a better resolved velocity profile is 
needed in order t o fit the exponent of the rotation velocity profile. 

Very recently, iTobin et al.l (|2012bl) reported the detection of a Keplerian disk 
around the Class protostar L1527 based on SMA and CARMA observations. 
This system is known to be edge-on with infrared images revealing on very small 
scales a 60 AU thick dark lane that suggests a disk seen in absorption against 
scattered light from the central protostar. The SMA continuum observations with 
0.25" angular resolution resolve a flattened structure coinciding with the infrared 
dark lane, with a radius of 90 AU and a mass of 0.007± 0.0007 Mq. The "CO 2-1 
CARMA observations show a clear velocity gradient along the flattened structure. 
The angular resolution is about 1" (140 AU), but the location of the emission 
peak in each velocity channel was determined with higher precision, which allowed 
the authors to derive a velocity curve consistent with v^ot oc r~'^'^. Assuming the 
disk is rotationally supported, they derive a mass of 0.19 ±0.04 Mq for the central 
protostar, and then an age of 3 x 10^ yr from an estimate of the mass accretion rate 
(6.6 X 10~^Mq yr"-*^) based on the measured accretion luminosity and assuming 
i?i = 1.7 Rq. This age is likely to be overestimated because it assumes a constant 
accretion rate. The dynamical timescale associated with the outflow is indeed 
an order of magnitude shorter. Like for NGC 1333 IRAS 4A, higher-angular- 
resolution molecular-line observations will be decisive to fit the exponent of the 
velocity profile and fully demonstrate the Keplerian nature of this disk. L1527 has 
the advantage of being a factor of ^ 2 closer. 



5.4 Class I protostellar disks 

As was seen in Sect. 15.31 the detection of rotationally supported disks in the Class 
phase is challenging because the envelope is still prominent at this stage and its 
structure needs to be well understood in order to disentangle the contribution of 
a disk to the detected emission. From this point of view, the Class I phase is a 
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priori more favorable because the mass of the central object (star plus disk) now 
dominates the system. 

The kinematics of the young (Tboi = 95 K) and luminous (Lboi = 33 Lq) Class I 
protostar L1551-IRS5 in Taurus was studied in C^ ^O 1-0 with the Nob eyama 
Millimeter Array (NMA) {HPBW 2.8" x 2.5") by iMomose et all ()l998t ). The 
continuum emission is elongated roughly perpendicular to the outflow and a clear 
velocity gradient is detected (see their Figs. 1 and 3). Its direction makes an 
intermediate angle to the outflow axis, which is interpr eted as resulting from a 
combination of infall and rotation. Momose et al. I ()l998h interprete the position- 
velocity diagrams in terms of a flattened envelope inclined to the line of sight. 
Infall produces a velocity gradient along the minor axis of the continuum emission 
and rotation along the major axis. They derive an infall velocity Vinf oc r~^'^. 
The rotation velocity is found to be Urot (x r~-^ in the outer parts (r > 5") and 
consistent with both and r""'^ in the inner parts. The authors estimate that 
contamination by the outflowing gas is negligible and derive a mass infall rate of 
6 X 10"^ Mq yr~^. They estimate the centrifugal radius to be ~ 160 AU, the 
radius were the infall and rotation velocities are approximately equal. A similar 
analysis was done for a small C^^O 1-0 survey of 2 Class and 6 Class I protostar s 
in Taurus performed with the NMA and OVRO (jOhashi et al.lll997l:l0hashilll999l ). 
The infall and rotation velocity profiles were derived from P-V diagrams along the 
minor and major axes, respectively. They found 3 sources with both infall and 
rotation signatures. From an extrapolation of the derived velocity profiles, they 
find centrifugal radii in the range 100-170 AU. These radii may characterize the 
rotationally-supported disks in the Class I phase but it should be kept in mind 
that they result from an extrapolation and are therefore relatively uncertain. It is 
in any case not a direct proof of the existence of such disks. 

The Class I sources included in the PROS AC sample (see Sect. 15. 3p were also 
observed with SMA in IICO+ 3-2 at a resolution of ~ 3", i.e. about 400 AU 
(jBrinch et al.ll2007t iLommen et al.ll2008l: I Jorgensen et al.ll2009l ). The HCO+ emis- 
sion was found to be elongated like the continuum emission for 4 out of 10 sources. 
For the other sources, it may be contaminated by the outflow. The P-V diagrams 
along the direction perpendicular to the outflow were compared to Keplerian ve- 
locity profiles for these 4 sources. The authors find such velocity profiles to be 
consistent with the data, suggesting that they may trace rotationally-supported 
disks (of sizes ^ 500 AU). However, these comparisons are not a proof that the 
"fit" is unique. Other velocity profiles may also be in reasonable agreement. There 
is a clear need for better resolved rotation velocity profiles for all Class I sources 
discussed up to here. 

HH 111 is a young (Tboi = 78 K) Class I system d riving a hi ghly collimated jet 
in L1617 in Orion (see Fig. 1 of Reipurth et al.lll999l ). iLed ( 201ll ) resolved the inner 
continuum emission as a flattened structure (deconvolved FWHM 240 x 120 AU) 
perpendicular to the jet. The envelope is rotating, a s traced in C^^O 2-1 emission 
with SMA [HPBW - 1.2", i.e. 500 AU. lLe3[2O10l) . The P-V diagram along the 
direction perpendicular to the jet shows a clear velocity pattern with v-^ot 
for R > 2000 AU and Vrot oc R-°-^ for R < 2000 AU (see their Figs. 4 and 
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5). This P-V diagram and the one along the outflow axis are fitted with an 
LTE radiative-transfer model, which provides a good agreement, apart from some 
contamination by the outflow in the latter (see their Fig. 7). The specific angular 
momentum is found to be uniform (8 fO~'^ km s~^ pc~^) over the range 2000- 
8000 AU, and decreasing toward the center at smaller radii. At higher- angular 
resolution with SMA {HPVW = 0.6", i.e. 240 AU), ^^CO 2-1 is, apart from 
some outflow contamination, well fitted by the rotation velocity profile der ived 
abov e from C^^O and suggests the presence of a resolved Keplerian disk (|Lee 
20111 ). From the continuum emission the author estimates a disk mass of 0.14 Mq 



and from the rotatio n ve l ocity profile (oc r~'^-^) in the inner part a stellar mass of 
1.3 Mq. In addition, Lee ( 2010[ ) argues that the compact emission probed with SO 
may trace an accretion shock at 400 AU. The picture emerging from these results is 
that HH 111 has an outer (2000-8000 AU) envelope collapsing with conservation of 
specific angular momentum, a transitional, sub-Keplerian, collapsing inner (400- 
2000 AU) envelope, and a rotationally-supported disk of radius 400 AU. One puzzle 
in this picture is that the rotation velocities are found to dominate over the infall 
velocities up to at least 8000 AU in the envelope. 

Finally, let us discuss the structure of L1551- NE, a young (Tboi — 91 K) Class I 
binary system in Taurus. iTakakuwa et al.l (|2012| ) observed it in 0.9 mm continuum 
emission and in C^^O 3-2 with SMA at very high angular resolution {HPBW ~ 
0.7", i.e. 100 AU). The continuum emission is ring-like, with a strong peak in 
the center, and is interpreted as consisting of (unresolved) circumstellar disks and 
a (resolved) circumbinary disk. The C^^O emission is modeled with both a thin 
Keplerian or infalling disk. The best-fit is obtained for a Keplerian-disk velocity 
structure without infall. From this modeling, the authors derive a circumbinary 
disk radius of 300 AU, a total (circumbinary -|-circumstellar) disk mass of 0.05 Mq 
and a stellar (binary) mass of 0.8 Mq. 

The presence of rotationally supported disks is proven rather convincingly for 
the last two sources, HH 111 and L1551 NE, with radii of 300-400 AU. Both 
sources have a bolometric temperature close to the threshold separating Class I 
from Class sources (70 K, but with the caveat that Tboi is sensitive to the 
inclination). The presence of extended Keplerian disks around such very young 
Class I sources strongly suggests that these extended disks were formed during 
the Class phase already. This holds true only if these sources are as young as 
suggested by their bolometric temperatures. The PROSAC sample contains good 
disk candidates, but higher-angular resolution observations of optically thin tracers 
are needed to confirm the Keplerian rotation in these objects. Overall, the studies 
mentioned in this section suggest Class I disk radii in the range 100-400 AU and 
a disk-to-star mass ratio in the range 1-10%, somewhat higher than for Class II 
sources (median ratio 0.9%, see Sect. 15. 1| ). 



5.5 Class II protoplanetary disks 



While Keplerian rotation in circumstellar disks is still laborious to probe in the 
Class and I phases, it is observationally well established around T Tauri stars. 
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This is for instance well shown by the CO 2-1 survey of T Tauri disks in Tau- 
rus carried out wi th the PdBI at an angular resolution of ~ 0.7" (100 AU) by 



bimon et alJ (I2OOOI) . They perfor med a minimizat i on of a parametric model o f 



a disk in hydrostatic equilibrium ( Dutrev et al. 1994 : Guilloteau fc Dutrevlll998l ). 
The velocity field is parametrized as Vrot oc r~°' and the best fit for all sources is 
obtained for a = 0.5, with a precision better than 10%. This allows for a precise 
measurement of the stellar masses. 

Three recent surveys of PMS circumstellar disks were made in continuum emis- 
sion: 14 low- and intermediate mass PMS stars (10 in Taurus) with CARMA 
{HPBW ^ 0.7", i.e. 100 AU, llsella et al.l l2009l) . 17 PMS stars in Ophiuchus 
with SMA {HPBW = 0.3", i.e. 40 AU. lAndrews et allliooi l2010h . and 23 PMS 
stars in Ta urus-Auriga with the P dBI at two frequencies {HPBW down to 0.4", 
i.e. 60 AU. iGuilloteau et~aL 2011). For the first two datasets, the surface density 
profile was derived by fitting a parametric model of a 2D, fiared density structure 
based on a similarity solution descri bing the viscous evolution of an accretion disk . 
with a viscosity scaling as oc i?^ (iLvnden-Bell fc Pringld l I1974 iHartmann et al 



1998[) . The radiative transfer was performed with a two-layer approximation for 



the first study and in full 2D for the second study. The dust opacity coefficient k 
and the gas-to-dust mass ratio were assumed to be uniform. The best-fit models 
of the first study yield disk radii in the range 90-320 AU, initial radii Ri ^ 25- 
40 AU (63% of the mass contained within 90% within 2i?i ) , and initial masses 
in the range 0.05-0.4 Mq. Making a number of assumptions, Isella et al. ( 20091 ) 
estimate the specific angular momenta required for parent dense cores to pro- 
duce disks with initial radii within the range they derive. They find jcorc ~ 0.8- 
4 10~* km s~^ pc, 10 times smaller than the observed specific angular momenta 
of N2H+ dense cores of mass 1-10 Mq. They conclude that only 10% of the core 
specific angular momentum was conserved during the collapse phase. Follow ing 
a similar procedure with slightly different assumptions, Andrews et al.l ( 2010[ ) de- 
rive j < 1.6-26 10"'' km s~^ pc for their Ophiuchus sample. They argue that 
this matches the lower range of N2H+/NH3 dense cores, and that the agreement 
may even be better if the observed dense-core j2D overestimate the true jsd (see 
Sect. 13. 6p . It is therefore unclear whether these disks measurements imply a loss 
of angular momentum during the collapse phase or not. 

The PdBI sample was fitted with two models, a truncat ed power-law model and 
a viscous disk model similar to the previous two studies. IGuilloteau et al.l (|2011l ) 
find that both models fit the data equally well. The measurements were done at 
two frequencies, allowing the dust opacity index /3 (k cx i/^) to be determined 
across the disks. /3 is found to increase from ~ at center to 1.7-2 at the edges, 
suggesting that grains are larger toward the center. The disk outer radii are found 
in the range 14-600 AU, with a hint of increase with stellar age, which would be 
consistent with a viscous evolution. Based on the viscous disk model, the authors 
derive initial disk radii smaller than 100 AU. 

All three studies find initial and current outer radii of PMS disks somewhat 
smaller than the disk radii derived for the Class I objects (~ 100-400 AU, see 
Sect. 15. 4p . This may suggest that the Class I "disks" are not yet fully rotationally 



28Angular momentum transport during the formation and early evolution of stars - EES2012 



supported disks or that they are pseudo-disks (see definition in Sect. 15. 3p . Alter- 
natively, inferring disk evolution from current surface density profiles of observed 
PMS disks based on a viscous disk model may not be correct. Finally, the samples 
are still small, especially for the Class I phase, and may not be free of biases. 



5.6 Summary 

The theoretical outcome of dense core collapse is uncertain in terms of disk forma- 
tion: the existence of disks and their properties depend strongly on the magnetic 
field configuration and on the role played by non-ideal MHD effects. The exis- 
tence of rotationally supported disks in the Class phase is not fully established 
observationally, even if two sources turn out to show good Keplerian-disk candi- 
dates. Higher angular resolution is needed to distinguish Keplerian disks from 
pseudo-disks or envelope small-scale structures. "Disks" associated with Class I 
protostars have typical radii of 100-400 AU and disk-to-star mass ratios of 1-10%. 
The presence of a rotation velocity field scaling as r~°-^ is well demonstrated in 
a few cases, possibly very young Class I objects, suggesting that these disks were 
already formed during the Class phase. The scaling as r"*^'^ may however not be 
an unambiguous sign of Keplerian rotation. Keplerian rotation is well established 
in disks around PMS stars. The small radii of these PMS disks question the Class I 
disks being already Keplerian. The latter may correspond to pseudo-disks. 



6 Rotation of jets 



6.1 Properties and origin of jets 

Jets and outflows are ubiquitous in star formation. Mos t YSOs wher e accretion or 
infall is occuring are associated with a jet or an outflow (|Cabritll2002h : Class and 
I protostars, which have strong or residual infall in their envelope plus possibly 
an accretion disk, and Class II PMS stars, which have accretion disks. A strong 
correlation between ejection and accretion was established in Cla ss II sources, sug 



gesting that the ejection process is driven by disk accretion (,Cabrit et al.lll99' 



The ratio is about 0.1-0.2 in Class II sources (|Cabritll2007h . Jets from Class II 



sources have narrow opening angles (a few degree) beyond 50 AU and the sources 
have no dense envelopes to confine them, which means that they must be intrin- 
sically coUimated. It is believed tha t this happens throu gh MHD self-collimation 



along rotating open field lines (e.g., iFerreira et al.l 120061 ). Jets are thought to be 



driven by magnetocentrifugal acceleration, which implies that they remove an- 
gula r momentum from the accretion disk and / or the stellar magnet osphere (see, 
e.g.. lK6nigl fc Pudritj[2000llShu et"aIll2000l:IZanni fc Ferreirall2013l and also the 
contribution of J. Ferreira in this volume). 

Three main magnetocentrifugal mechanisms have been proposed fo r steady jets 
(see, e.g.. iKonigl fc Pudritz|[200ol: IShu et "alll2000l: iFerreira et al.lbood for details): 
the stellar wind model which launches the jet close to the stellar surface, the 
X-wind model for which the jet is launched from a narrow region close to the 
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radius Rx of the disk where it is truncated by the stellar magnetosphere, and the 
extended disk wind model for which the launching occurs over a wider range of 
disk radii (> Rx)- Constraints on the launching mechanism, such as the radius 
where it occurs and the magnetic lever arm braking the rotating disk/star, can 
be derived from the location of the jet in a diagram displaying the local specific 
angular momentum of the jet {Rv^{R), with Vd,(R) the toro idal velocity) versus 
its poloidal velocity (vp{R)) (see Fig. 2 of lFerreira et al. 2006 ). This motivates the 
sea rch for rotation s ignat ures in jets. 

Anderson et"al] (|2003l ) computed a practical form of the equation relating the 



launching radius Rq to the poloidal and toroidal velocities of the jet, Vp{R) and 
v^{R), at a distance R from its axis and far from the launching region (their 
Equation 5): 



Ro ~ 0.7 AU 



R 



10 AU 



2/3 



V4,iR) 

10 km s" 



2/3 



Vp{R) 



100 km s" 



-4/3 



1 Mp 



1/3 



(6.1) 

As stated by lAnderson et al.l ( 20031) . this equation is an approximation and is 
valid only when the kinetic energy of the jet is much greater than the gravitational 
binding energy at the launching region. It is not valid anymore when the magnetic 
lever arm is small, in which case more terms have to be taken into account (see 
th eir Equation 4 ) . Th is is the reason why the curves of constant Rq in Fig. 2 
of lFerreira et al.l ( 2006f) are curved toward low values of the magnetic lever arm 
(A) while Equation 16.11 would produce straight lines. After using Equation 16. 1[ 
it should thus always be verified a posteriori that the gravitational energy at the 
launching region is negligible compared to the kinetic energy of the jet, otherwise 
the complete equation has to be used. 



6.2 Jet rotation signature 

In the absence of rotation, spectra taken at symmetric position on each side of the 
axis of an axisymmetric jet should have the same centroid vel ocity. If the jet is ro - 
tating, there should be a velocity shift between both spectra. iPesenti et al. I (|2004l) 



present predictions for an optical line ([O I] at A6300) based on a self-similar MHD 
disk-wind solution and taking into account projection effects (inclination angle i) 
and a finite angular resolution. They investigate two ionization profiles, Xe (x ^ 
and Xc on Rq, Rq being the radius in the disk. Let's call vr and V-r the centroid 
velocities of spectra taken at symmetric positions at distances +/—R on each side 
of the jet axis, Vp{R) and v^{R) the poloidal and azimuthal velocities of the jet, 
and VsMtt (R) the velocity offset between both spectra. Under the assumption that 
only one streamline dominates the emission along each line of sight, we have the 
following equations: 

^p(^) = 7— ^ ' (6-2) 



(6.3) 
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Equation 16.21 assum es that fn ~ "^ ^zi i-e . that the radial velocity Vr « Vp. As 
shown in Fig. 3b of iPesenti et al.l (|2004l ). v^{R) can be underestimated if the jet 
section is not well resolved (for R < 2i?{,eam = HPBW). This implies that only 
the external radius of the launching region can often be reliably determined. v^{R) 
is also sensitive to the distribution of the line emission across the jet axis because 
in reality several streamlines contribute along each line of sight. If the outer, 
slower streamlines dominate the emission, then (i?) will also be underestimated 
(see the difference between the predictions for oc and Xe oc i?o in Fig. 3b of 
Pesenti et al.ll20()3 ). As a result, a careful choise of the tracer and a high angular 



resolution (better than 5 AU, i.e. 0.04" at 125 pc) are needed to fully probe the 
rotation profile of the jet and test MHD launching models. 



6.3 Jet rotation in the Class II and I phases 



Indications of jet rotation based on measurements of transverse velocity shifts 
were found in a few Class I and II sources (see Sec t. I6.6p. One of t h e first claimed 
evidences for rotation in a j et was pr esented in Bacciotti et al.l ( 20021) for the 
Class II object DC Tan (see iDavis et al. 2000 for an earlier claim in HH 212). 
The authors observed with a slit parallel to the jet, at seven offsets across the jet, 
in optical forbidden lines ([O I], [N II], and [S II]) with an angular resolution of 
~ 0.1", i.e. about 14 AU. They applied two methods, multiple Gaussian fitting 
and cross-correlation of spectra at symmetric positions with respect to the jet axis, 
both with an accuracy of 5 km (see their Fig. 1). They measure velocity 
shifts Wshift ^ 5-20 ±5 km s^^ across the jet (see their Fig. 2) and conclude that the 
jet is rotating with ^ 6-15 km s^^ for distances from the central star between 
0.075" and 0.4", i.e. 11 and 56 AU (in projection). In the framework of the jet 
launching paradigm described in Sect. 16. 1[ the r adius of the launch ing region is 
i?o ^ 3 AU for the external par t of the atomic j e t IPesenti et aL 2004 , with a mor e 



sophisticated analysi s than i n iBacciotti et al.l |2002| and I Anderson et al.l 120031 ) . 
In addition, Baccio tti et al.l (j2002n estimate that the flux of angular momentum 
carried by the jet represents about 60% of the angular momentum loss rate needed 
in the disk at Rq for accretion to occur at the observed rate. Further observations 
with HST/STIS in t he near-UV c onfirmed the sign and magnitude of the velocity 
gradient in DC Tau ( Coffev et al.ll2007[ ). Signatures of rotation were also found for 
the molecular (H2) jet of the Class I protos tar HH 26 with ob servations performed 
in H2 1-0 S(l) with VLT/ISAAC (Chrvsost'o^ou et al.ll2008h . The authors derive 
i?o ~ 2-4 AU and estimate that the transported angular momentum flux represents 
about 70% of the loss rate needed for accretion. These two examples suggest that 
jets significantly contribute to the removal of angular momentum from accretion 
disks at the launching radius, both in the Class I and II phases. 

Altough the previous examples suggest that jet rotation is observationally well 
established, contradictions found in other sources raise the question whether these 
tiny velocity shifts really trace rotation. Using optical lines with HST/STIS, 



Woitas et al l (|2005h reported the detection of rotation in the jet of RW Aur based 
on such velocity shifts, but with a large dispersion (see their Figs. 4 and 5). How- 
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ever, the disk of RW Aur was found to rotate in the oppo site direction, based 



on CO 2-1 observations with the PdBI ( Cabrit et al. 20061) . which casts serious 



doubts about the interpretation of the velocity shifts in the jet as tracing rota- 



tion. Further observations of the jet were performed in the near-UV (jCoffev et al 
2OI2I) . Velocity shifts in the approaching jet lobe were found to be consistent with 



the disk rotation, thus in contradiction with the velocity shifts measured in the 
optical. But no velocity shift was apparent in the receding jet, and no velocity 
shift was found again in the approaching jet six months later. The detection of jet 
rotation in RW Aur has therefore not been secured yet. 



6.4 Jet rotation in the Class phase 

Evidence for jet rotation has been claimed for several Class protostars. For 
instance, velocity gradients across the jet of NGC 1333 IRAS 4A2 were reporte d 



in SiO 1-0 with the VLA {HPBW ~ 1.5-2", i.e. 350-470 AU. IChoi et al.ll201lD . 
The velocity gradient has the same direction for both lobes and is seen along 8 
cuts perpendic ular to the jet. I n addition, the direction of rotation is the same as 
for the "disk" ( Choi et al. 2010l see Sect. 15. 3p . One caveat concerning the analysis 
is that the velocity gradient is inferred from simply connecting "blobs" in the P-V 
diagrams and is therefore somewhat subjective (see Fig. 2 of Choi et al. 201 ih . 



The proper motion of the H2 outflow was measured by I Choi et aL ( 20061) who 



derived a poloidal velocity ^ 70 km s~ . It was assumed to be the same for 
the SiO component. The radius of the SiO jet is found to increase linearly and 
the angular s peed of the oute r layer to decrease with distance from the protostar 
(see Fig. 3 of Choi et al. 201 1[) . The specific angular momentum slightly increases 
with distance, which led the authors to suggest that the angular momentum in- 
jection at the basis of the jet decreases with time. At least, there seems to be no 
significant loss through interaction with the ambient medium. In the framework 
of the magnetocentrifugal launching paradigm, the derived poloidal and azimuthal 
velocities imply a launching radius Rq ^ 2 AU (see their Fig. 5), which favors the 
extended disk-wind model over the X-wind and stellar-wind models, and a mass 



ejection efficiency as traced by the SiO component 



Mac 



1%. The latter 



is only a crude estimate and needs to be verified with a measurement of the mass 
ejection rate. 



6.5 Outflow rotation 

So far, we discussed claims for rotation in jets. Rotation has also been searched for 
in molecular outflows. Molecular outflows are less coUimated than jets which typ- 
ically have a full opening angle < 10°. The former may consist of cloud material 
entrained by the latter, or material launched at larger radii than je ts and maybe 



appe a ring at the stage of the first hy drostatic core already (see, e.g. Ma chida et al 
l2008t iHennebelle fc Fromand 120081) . The molecular outflow of CB26, a Class I 



object in Taurus-Auriga, was mapped with the PdBI in CO 2 -1 with an angu- 
lar resolution of ~ 1.5" (about 210 AU. iLaunhardt et aLllioOol ). It is small and 
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well collimated, and it has the same direction as the HH objects associated with 
CB26. A clear velocity gradient is detected orthogonal to the outflow axis (see 
their Fig. 4). The gradient h as the same orientation as the one measured in the 
disk. Launhardt et al. I (|2009l) perform an empirical, parametric modeling of the 
disk and outflow interferometric data, also using constraints from the spectral en- 
ergy distribution and near-infrared maps. They find that a rotating outflow is 
consistent with the observations (see their Figs. 5 and 6), even if the exact shape 
of the rotation velocity profile across the outflow is not well constrained. They 
derive a mass and total angular momentum of the outflow two orders of magni- 
tude smaller than the mass and angular momentum of the disk, and a specific 
angular momentum in the outflow similar to the one in the disk. They conclude 
that it will take about 1 Myr for the outflow to dissipate the mass and angular 
momentum of the disk, which is comparable to (although somewhat shorter than) 
the statistical dispersion timescale of disks in the Class II phase. In the frame- 
work of the magnetocentrifugal paradigm, the measurements imply a la unching 
radiu s Rq ~ 5 ± 4 AU, the uncertainty being related to the inclination (jCabrit 
2009f) . Although the in terpretation in terms of rotating outflow seems plausible. 



Launhardt et al. do not exclude alternative interpretations, such as the 



presence of two unresolved outflows or jet precession. 



6.6 Statistics 



If the interpretation of the velocity gradients across the jets in terms of rotation is 
correct and if the magnetocentrifugal paradigm is the right one, then the current 
m easurements are all consistent wit h the extende d disk-wind model (see Fig. 3 
of Ferreira et al. 2006 and Fig. 5 of Cabrit 20091 and the discussions therein). 
The ratios of mass ejection to mass accretion rates inferred from the analysis are 
■j^sj- < 0.1-0.3 ( Ferreira et al. 20061 ). They are consistent with the ratios directly 



measured for Class II sources (- 0.1-0.2. ICabr"i^l2007l) . 



A summary of most measurement s of rotation in jets or outflows is presented 
in Table [2] (adapted from Table 1 of iBacciottil l2009l ) . Out of 13 Class 0, I, and 
II sources which were observed, 10 show a velocity gradient perpendicular to the 
jet/outflow axis. For 5 of them, the orientation of the velocity gradient is consistent 
with the disk rotation, which makes the interpretation in terms of jet rotation 
very plausible. However, in 4 cases, there is a disagreement either between the 
jet/outflow and disk directions of rotation (RW Aur, HH 212), or between the 
two jet/outflow lobes (HH 212), or no velocity gradient is found while the high 
inclination of the system to the line of sight {i ^ 90°) should have been favorable for 
a detection (HH 30, HH 111). These cases cast some doubts about the rotation 
interpretation. The detection of jet/outflow rotation is therefore not yet fully 
established. 
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Table 2. Measurements of rotation in jets/outflows. 



Jet/ 


YbO 


A 


Obs. 


AV 


Bipolar 


Jet/ 


Ref.*^ 


outnow 


ciass 


range 


mode 




lobes 


disk 




DG Tau 


II 


VIS 


HST II 


Y 


- 


agreem. 


1, 2 






VIS 


HST _L 


Y 


- 


agreem. 


3 






NUV 


HST _L 


Y 


- 


agreem. 


3 


RW Aur 


II 


VIS 


HST II 


Y 


agreem. 


disagr. 


4, 5 






VIS 


HST _L 


Y 


agreem. 


disagr. 


6 






NUV 


HST _L 


Y/N 


? 


agreem. 


7 


TH 28 


II 


VIS 


HST _L 


Y 


agreem. 


- 


6 






NUV 


HST _L 


Y 


agreem. 


- 


3 


CW Tau 


II 


VIS 


HST _L 


Y 


- 


agreem. 


3 


HH 30 


II 


VIS 


HST _L 


? 




inconcl. 


6, 8 






mm 


PdBI 


N 




inconcl. 


8 


CB26 


I 


mm(CO) 


PdBI 


Y 


agreem. 


agreem. 


9 


HH 26 


I 


NIR 


VLT ± 


Y 






10 


HH 72 


I 


NIR 


VLT ± 


Y 






10 


HH 111 


I 


NIR 


GEMINI 1 


N 




disagr. 


11 


HH 34 


I 


NIR 


GEMINI ± 


U 






11 


HH 212 





NIR 


UKIRT II 


Y 


disagr. 


agreem. 


12, 13 






NIR 


GEMINI ± 


Y 


compat. 


agreem. 


11 






mm(SiO) 


PdBI 


Y 


disagr. 


disagr. 


14 






mm(SiO) 


SMA 


Y 


agreem. 


agreem. 


15 


HH 211 





mm(SiO) 


SMA 


Y 


agreem. 


agreem. 


16 


IRAS 4A2 





mm(SiO) 


VLA 


Y 


agreem. 


agreem. 


17 



Notes: adapted from Table 1 of lBacciottil (|2009[ ). ° Observing mode. "||" and "_L" 
mean observations with a slit parallel or perpendicular to the jet axis. ^ "Y"/"N" 
mean that a velocity gradient is/is not detected. "U" means that the emission 
was unresolved. "agreem." /"disagr." means that a velocity gradient is found in 
both lobes with the same/opposite orientation. A dash means that only one lobe 
was observed. "agreem." /"disagr." means that the jet/outflow and the disk are 
found to rotate in the same/opposite direction. A dash means that the direction 



of the disk rotation is not known. References: 1: Bacciotti et aL (2002), 2 



12: 



Davis et al 



Testi et all (l2002h. 3: [Coffev et all (|2007h. 4: IWoitas et al. (20051 5: ICabrit et al 



(l2006l). 6: [Cofltev et all (l 2004h. 7: Icoffev et al.l (l2012l) 8: [Petv et all (l2006l). 9 



Launhardt et all (l2009l) . 10: IChrvsostomou et all (l2008l) 11: ICoffev et al.l (l201ll) 



Lee et all (|2o7)8h 16: iLee et al.l (|2007h . 17: Ichoi et aD (|201lh 



(l2000l). 13: IWiseman et al.l (l200lh. 14 : ICodella et alJ (l2007l). 15: 
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6.7 Alternative interpretations 

Several alternative interpretations instead of rotation have been proposed to ex- 
plain the velocity gradients across jets/outflows. Beyond the Alfven surface, a few 
AU above the disk surface, a jet may undergo various kind of instabilities which 
can produce differential velocities on the order of the sound speed 10 km s~^), 
similar to the v elocity shifts measured across the jets and interpreted as rotation 
( Bacciottilliooi ). ] t is however likely that such instabilities would produce random 



or periodic velocity gradients as a function of distance from the center, and not a 
regular pattern as is observed in a few sources, with the same orientation in both 
lobes (e.g. TH 28, CB26, NGC1333 IRAS 4A2). 



The effect of jet precession was investigated numerically by iCergueira et al 



(200i). They performed 3D simulations of a non-magnetic jet with variabilities in 



ejection direction (precession) and jet velocity (intermittence) , and with or with- 
out rotation. An atomic and ionic network was included in order to compute 
line profiles. They found that only precession (with intermittence) or only rota- 
tion produce very similar velocity shifts, implying that such velocity shifts do not 
necessarily trace rotation. 

Soked (l2005l) proposed a phenomenological model of jet interaction with a 



warped disk to explain the velocity shifts seen across jets. In this model, the 
jet interacts with the ambient gas above the disk surface and is more slowed down 
on the warped side of the disk than on the opposite side. In this way, one side 
of the jet is slower than the other, which produces a velocity shift between both 
sides, not due to rotation. In this model, the velocity shifts would not be seen for 
inclinations close to 90°, i.e. a jet propagation in the plane of the sky. This would 
explain why no velocity shift is seen in the jets of HH 30 or HH 111. 

Finally, if the ambient medium is not uniformly distributed, its interaction 
with the flow could create asymmetries in the propagation of the jet, which woul d 
in turn produce velocity shifts between both sides of the jet ( Bacciottl 2009f ). 



However, persistent velocity patterns with increasing distance from the star could 
a priori not be produced in this way. 



6.8 Summary 

The following requirements are needed to convincingly detect rotation in jets or 
outflows. First of all, the kinematic signature, a velocity gradient perpendicular 
to the jet axis, has to be consistent all along the jet length. Second, the signature 
should be consistent between both lobes and with the disk rotation. Finally, the 
velocity profile should be spectrally resolved. In light of this, the detection of jet 
rotation is statistically not yet secure. Doubts remain because some cases present 
inconsistencies between the putative directions of rotation of the disk and the jet. 
In addition, the absence of velocity shifts in nearly edge-on systems such as HH 30 
or HH 111 is puzzling since they should be the most favorable systems for the 
detection of jet rotation. Alternative interpretations ~ jet instabilities, precession 
coupled with modulation of jet ejection, interaction of the jet with a warped disk. 
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asymmetric shocking - let open the question whether the measured velocity shifts 
really probe jet rotation or not. 

Nevertheless, if the interpretation in terms of rotation is valid, then the current 
constraints on the poloidal and azimuthal velocities are consistent with extended 
disk- wind models with an external launching radius of a few AU. The Class II jets 
would seem to account for a large fraction (60-70%) of the loss rate of angular 
momentum in disk that is required at the launching radius to enable accretion at 
that radius. 

7 Conclusions 

7.1 Reliability of rotation signatures 

The previous sections introduced the methods used to probe rotation in the dif- 
ferent phases of star formation. On large scales 1 pc), the velocity gradients 
measured in molecular clouds imply a specific angular momentum scaling as _R~^'^ 
if they are interpreted in terms of rotation, possibly meaning that the contraction 
of molecular clouds occurs with a loss of angular momentum. However, this rela- 
tion may simply be related to the scaling properties of turbulence, suggesting that 
the velocity gradients trace turbulence vorticity rather than pure, well-ordered 
rotation. On the smaller scales of dense cores 0.1 pc), the correlation is simi- 
lar, but solid-body rotation seems to be rare, and the velocity structure of dense 
cores in protoclustcrs is relatively complex. At a later stage, the velocity gradi- 
ents found in protostellar envelopes are often not orthogonal to the outflow axis 
and the rotation signature may be ambiguous if the envelopes are prolate rather 
than oblate or spherical. However, there are a few "well-behaved" cases for which 
the signature of rotation is relatively robust. On small scales (< 100 AU), the 
existence of rotationally-supported disks in the Class phase is observationally 
not established yet. In the Class I phase, rotationally-supported disks have been 
claimed to be detected but there is a size issue compared to the Class II disks and 
the former could be "pseudo-disks" . The most robust detections of rotation have 
been made for the Class II phase, where the kinematic structure of the detected 
compact structures leaves little doubt that they are rotationally-supported disks, 
i.e. structures in Keplerian rotation. The ubiquitous jets found at all stages of 
star formation from the Class to the Class II phase are expected to rotate if 
the current paradigm of magnetocentrifugal launching holds. Several detections 
of jet rotation have been claimed, but they are statistically not secure yet. If the 
interpretation of the measured velocity gradients in terms of rotation is valid, then 
these detections imply a high efficiency of angular momentum removal from the 
disk by the jet at a distance of a few AU from the protostar. 

7.2 Evolution of angular momentum 

As a summary of this review. Fig. [8] shows the distribution of specific angular mo- 
mentum as a function of radius for different evolutionary stages - from molecular 
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Fig. 8. Specific angular momentum as a function of radius for different evolutionary 
stages or different scales of a star forming region. The specific angular momentum is the 
local value (i.e. the product of the rotation velocity and the radius) for all categories 
except for the binaries, the solar system, and the Sun for which it truly corresponds to the 
mean angular momentum per unit of mass. For most categories, the radius is the radius 
at which the rotation velocity was measured or fitted and does not necessarily correspond 
to the outer radius of the object. The two solid lines are intended to guide the eye and are 
not least-square fits to the data. The vertical dashed lines mark the approximate position 
of breaks in the distribution of angular momentum and the horizontal one indicates the 
typical specific angular momentum during the protostellar collapse phase. References are 
listed in Table El 



clouds to the Sun - or different scales of a star forming region - e.g., protostellar 
envelope, circumstellar disk, jet. Referenc es are lis t ed in Table |3l This f i gure i s 
par tly based on diagra ms published by Oha shi et al.l (jl997l ). lBelloche et al.l (|2002l ). 
and lChen et aP (|2007h . There are several caveats to this plot. First of all, for most 
categories of objects, the angular momentum was derived assuming that velocity 
gradients do trace rotation, which, as was seen in the previous sections, may not 
always be true. Second, a correction for inclination was applied for most categories 
but not all (see Table [3]). Third, the plot is a mixture of an evolutionary diagram 
- where the radius can be viewed as probing the contraction state of an object - 
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Table 3. References for angular momenta shown in Fig. [S] 



Category 


Symbol 


,/ a 
M 


Corr.'' 


Rcf."^ 


Molecular clouds and dense cores 








Molecular clouds 


open square 


local 


no 


1 


i > J. J.3 U-CJ-ioC ^jVJ-L CO 




local 


no 


2 


N2H+ dense cores 


open triangle (down) 


local 


no 


3 


Ophiuchus cloud / cores 


filled triangle 


local 


no 


4 


Protostellar envelopes 










IRAM 04191 env. (CI. 0) 


purple star symbol 


local 


yes 


5-6 


Class envelopes 


red star symbol 


local 


no 


7 


Class and I envelopes 


black star symbol 


1 1 

local 


yes 


8 


HH 111 env. (CI. I) 


blue star symbol 


local 


ves 


9 


Circumstellar disks 










Class "disks" 


red circle 


local 


yes 


10-11 


Class I "disks" 


black circle 


local 


yes 


12-15 


Class II disks 


blue circle 


local 


yes 


16 


Jets 










Class II jet 


blue arrow 


local 


yes 


17 


Young binaries 










Class I binaries 


black plus symbol 


mean 


no 


7, 18 


T Tauri binaries 


blue plus symbol 


mean 


no 


7, 18 


Solar system 










Jupiter 


\ 


local 


yes 


19 


Solar system 


black 


mean 


yes 


19 


Sun 


blue 


mean 


yes 


20 



X iiiot) or mean [jj) specific angular momentum. ^ 
corre ction for inclination was a pplie d. ^ References: 1 : 
Goldsmith fc Arauillal (Il985l). 2: ICoodman et all (Il993h. 3: ICaselli et al.l (l2002l) . 



Notes: ° Local (wrot 
Flag indicating if the 



11: 



4: iBellochd (l2002l) . 5 
Chen et al.l (l2007l) 



Tobin et al 



J0rgensen e t aL. 

Bacciotti et 
Pinto et al.r(|201l[) . 




Belloche e t al' ('20021) . 6: iBelloche fc Andrei (120043^7 



Ohashi et al, (1997) 9: Ee3 (l2010l). 10: IChoi et al.l (I2OIOI) 



Brinch et all (120071). 13: iLommen et al 



Takakuwa et al.l (|2012h . 16: ISimon _e t al 
Chen (2013, priv. comm.), 19: [Alien. (1975 



and a snapshot of the structure of several objects - for which the diagram displays 
the angular momentum profile at a given time. Finally, the sample is certainly 
not complete and, in particular, upper limits for objects where no rotation was 
detected - in some cases maybe because of an unfavorable inclination, though - 
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are not shown. 

Despite these caveats, the overall shape of the distribution of specific angular 
momentum as a function of radius is likely real. Three regimes can be distin- 
guished. On large scales, the angular momentum of molecular clouds and dense 
cores appears to correlate well with the radius. The specific angular momentum 
is reduced by nearly two orders of magnitude over two orders of magnitude in ra- 
dius. On the one hand, as summarized in Sect. 13.71 this loss of angular momentum 
could be due to magnetic braking or gravitational torques, both transfering an- 
gular momentum to the ambient medium, or fragmentation - part of the angular 
momentum of the fragmenting structure being transfered into the orbital motion 
of its fragments. On the other hand, the correlation between "angular momen- 
tum" (derived from velocity gradients) and radius may not be directly related to 
rotation but rather to the properties of interstellar turbulence and the velocity 
gradients may trace the behaviour of turbulence vorticity. More details about the 
physical processes at work during this phase can be found in the introductory 
chapter of Hennebelle et al. in this volume. 

The second regime concerns the scales between ~ 100 AU and 0.03 pc (1.5 
orders of magnitude), which correspond to the scales over which a protostellar 
envelope is dynamically collapsing. In this regime, as already summarized in 
Sect. 14.61 the s pecific angula r momentum seems to be relatively constant, as was 
first noticed bv lOhashi et al ](IL997). This suggests that processes such as magnetic 
braking or fragmentation are relatively inefficient in reducing the specific angular 
momentum over this range of scales during the protostellar collapse. Note that 



Yen et al.l (j201ir) discuss a possible evolution in time of the angular momentum 



plateau in this regime. 

Third, the few datapoints shown below ~ 100 AU in Fig. [8] as well as the 
resolved structure of Class II disks (see Sect. 15.51) suggest that Keplerian rota- 
tion dominates the kinematics on these small scales, which implies that efficient 
mechanisms removing a large fraction of the angular momentum from the inner 
parts are at work to allow accretion to proceed. Processes contributing to the out- 
ward transfer of angular momentum within circumstellar disks are discussed by S. 
Fromang in this volume. In addition, the tentative detection of rotation in jets, 
if confirmed, indicates that they do extract from circumstellar disks a significant 
fraction of the angular momentum in the (small) launching region (see also the 
review of J. Ferreira in this volume). Fragmentation and the formation of binaries 
on these small scales can also significantly help solving the angular momentum 
problem by storing a large fraction of the angular momentum into the orbital mo- 
tions of the stars. Systems forming single stars, like our solar system, do however 
not benefit from this mechanism. 

Finally, if we extrapolate the average level of specific angular momentum at 
100 AU in circumstellar disks down to the radius of the Sun with a power-law 
consistent with Keplerian rotation (wiot(^) x i? cx i?°-^), then we obtain a value 
still three orders of magnitude higher than the specific angular momentum of 
the Sun. Processes that can reduce the specific angular momentum of a young 
star during its evolution (star-disk interactions, stellar winds) are addressed by J. 
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Bouvier and J. Ferreira in this volume. 



7.3 Will ALMA help us? 

ALMA is approaching its completion. The early results obtained in Cycle with 
about one third of the final number of antennae are already impressive and very 
promising for the future operation with the full array. With its high angular and 
spectral resolutions, ALMA will be a prime instrument to study the kinematics of 
star forming regions at very small scales with molecular line observations. 



According to ISchievenI (j2012l ). the finest angular resolution obtained at the 
highest frequency 900 GHz) will be about 0.005" with the full array in its most 
extended configuration, i.e. 0.7 AU at 140 pc. The spectral line sensitivity for 
these extreme conditions - a > 9000 K in 1 min at a spectral resolution of 0.3 
km/s, or > 400 K in 10 h - will not be sufficient to detect thermalized (non-masing) 
transitions, but lower-frequency observations are still promising: for instance, a 
sensitivity of 330 K for a spectral resolution of 2 km s~^ and an angular resolution 
of 0.03" - 4 AU at 140 pc - should be achievable at 150 GHz in 1 min, which 
translates into 13 K in 10 h. The capabilities of ALMA should thus allow to verify 
and firmly establish the presence of rotation in jets as well as fully resolve rotation 
profiles of disks around Class and I sources and thereby significantly contribute 
to our understanding of the formation and evolution of circumstellar disks. 
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